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About This Document 


This document contains specifications for a math library that takes advantage of the Single Instruction, Multiple Data 
(SIMD) instructions provided by the PowerPC Processor Unit (PPU) and the Synergistic Processor Unit (SPU) 
hardware of the Cell Broadband Engine”. By computing multiple results at one time, SIMD math functions allow 
programmers to obtain much higher performance from their PPU and SPU programs than would be possible from a 
corresponding traditional scalar math library. 


Audience 


This document is intended for system and application programmers who are interested in writing high-performance 
programs for the Cell Broadband Engine”. 


Version History 


This section describes significant changes made to each version of this document. 





Version Number & Date Changes 








v. 1.0 Created the initial document. 
November 6, 2006 


Related Documentation 


The following table provides a list of references and supporting materials for this document: 





Document Title Version Date 








C/C++ Language Extensions for the Cell Broadband 2.3 December 2006 
Engine Architecture 


ISO/IEC Standard 9899:1999 (C Standard) 


IEC Standard 60559:1989 (Standard for Binary Floating- 
Point Arithmetic) 


Synergistic Processor Unit Instruction Set Architecture 1.11 October 2006 


PowerPC Microprocessor Family: Vector/SIMD Multimedia 2.06c 
Extension Technology Programming Environments Manual, 
http://www-306.ibm.com/chips/techlib/techlib.nsf/techdocs/ 
C40E4C6133B31EE8872570B500791108 


PowerPC Architecture Book, http://www-128.ibm.com/ 2.02 
developerworks/eserver/library/es-archguide-v2.html 


Document Structure 


This document contains two chapters. The first is a SIMD math library overview, and the second is a specification 
describing the particular math functions that compose this library. 


Bit Notation 


Bit Notation 


Standard bit notation is used throughout this document. Bits and bytes are numbered in ascending order from left to 
right. Thus, for a 4-byte word, bit O is the most significant bit and bit 31 is the least significant bit, as shown in the 
following figure: 
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+ MSB 





MSB = Most significant bit 
LSB = Least significant bit 


Notation for bit encoding is as follows: 


e Hexadecimal values are preceded by 0x. For example: 0x0A00. 


e Binary values in sentences appear in single quotation marks. For example: ‘1010’. 


Byte Ordering and Element Numbering 


As shown in Figure 1, byte ordering and element/slot numbering is always displayed in big endian order. 


Figure 1: Big-Endian Byte/Element Ordering for Vector Types 









Byte 0 
(MSB) 


Byte 1 | Byte 2 [Byte 3 | Byte 4 | Byte 5 | Byte 6 | Byte 7 | Byte 8 | Byte 9 |Byte 10 [Byte 11|Byte 12|Byte 13 [Byte 14 [Byte 15 



































doubleword O doubleword 1 








word O word 1 word 2 


halfword 1 halfword 2 halfword 3 halfword 4 halfword 5 halfword 6 halfword 7 


char 8 | char 9 | char 10 |char 11 | char 12 | char 13 | char 14 | char 15 







halfword O 





charO | char 1 | char2 | char3| char4 | char5 | char 6 | char 7 

















Typographic Conventions 


In addition to bit notation, the following typographic conventions are used throughout this document: 





Convention Meaning 








courier Indicates programming code, processing instructions, register names, data 
types, events, file names, and other literals. Also indicates function and 
macro names. This convention is used only where it facilitates 
comprehension, especially in narrative descriptions. 


courier + italics Indicates arguments, parameters and variables, including variables of type 
const. This convention is used only where it facilitates comprehension, 
especially in narrative descriptions. 


italics (without courier) Indicates emphasis. Except when hyperlinked, book references are in italics. 
When a term is first defined, it is likely to be in italics. 


blue Indicates a hyperlink (color printers or online only). 
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1. Overview of the SIMD Math Library 


The PPU and SPU instruction sets include Single Instruction, Multiple Data (SIMD) instructions, which are similar to 
normal instructions but operate on more than one input simultaneously. Traditional math functions operate on a single 
input and are unable to take advantage of the speed and power of SIMD instructions. The SIMD Math Library 
contains SIMD versions of the scalar math functions described in the C99 standard, or /SO/IEC Standard 9899:1999 
(C Standard). This chapter provides specifications for these special PPU and SPU SIMD libraries. 


A. Library and Header Files 


The name of the SIMD library will contain the string simdmath. For example. on GNU/Linux the library will be called 
libsimdmath.a, or libsimdmath. so (for the shared library version). The simdmath.h system header file will 
contain type declarations and prototypes for the SIMD math functions. 


B. Functions Overview 


The functions that comprise the PPU and SPU SIMD math libraries are listed in Table 1. The functions that are listed 
as “non-standard” have no C99 counterpart. 


Names of the SIMD math functions are differentiated from their scalar counterparts by a vector type suffix appended 
to the standard scalar function name. For example, the SIMD version of £abs£ () , which acts on a vector float, is 
called fabsf4 (). Similarly, a SIMD version of a standard scalar function that acts on a vector double will have 42 
appended to the name. 


Table 1: SIMD Math Functions 





Function Name C99 Name Function Category Precision SPU/PPU 
absi4 abs integer word Both 
acosd2 acos trig double SPU 
acosf4 acosf trig single Both 
acoshd2 acosh trig double SPU 
acoshf4 acoshf trig single Both 
asind2 asin trig double SPU 
asinf4 asinf trig single Both 
asinhd2 asinh trig double SPU 
asinhf4 asinhf trig single Both 
atan2d2 atan2 trig double SPU 
atan2f4 atan2f trig single Both 
atand2 atan trig double SPU 
atanf4 atanf trig single Both 
atanhd2 atanh trig double SPU 
atanhf4 atanhf trig single Both 
cbrtd2 cbrt power double SPU 
cbrtf4 cbrtf power single Both 
ceild2 ceil rounding double SPU 
ceilf4 ceilf rounding single Both 
copysignd2 copysign other double SPU 
copysignf4 copysignf other single Both 
cosd2 COS trig double SPU 


SIMD Math Library Specification for Cell Broadband Engine™ Architecture, Version 1.0 


Overview of the SIMD Math Library 


SONY 


SONY 








Function Name 


cosf4 
coshd2 
coshf4 
divd2 
divf4 
divi4 
divu4 
erfcd2 
erfcf4 
erfd2 
erff4 
exp2d2 
exp2f4 
expd2 
expf4 
expm1d2 
expm1f4 
fabsd2 
fabsf4 
fdimd2 
fdimf4 
floord2 
floorf4 
fmad2 
fmaf4 
fmaxd2 
fmaxf4 
fmind2 
fminf4 
fmodd2 
fmodf4 
fpclassifyd2 
fpclassifyf4 
frexpd2 
frexpf4 
hypotd2 
hypotf4 
ilogbd2 
ilogbf4 
irintf4 
iroundf4 
isOdenormf4 
isOdenormd2 
isequald2 
isequalf4 


C99 Name 


cosf 

cosh 

coshf 
non-standard 
non-standard 
div 
non-standard 
erfc 

erfcf 

erf 

erff 

exp2 

exp2f 

exp 

expf 

expm1 
expm1f 

fabs 

fabsf 

fdim 

fdimf 

floor 

floorf 

fma 

fmaf 

fmax 

fmaxf 

fmin 

fminf 

fmod 

fmodf 
fpclassify 
fpclassify 
frexp 

frexpf 

hypot 

hypotf 

ilogb 

ilogbf 
non-standard 
non-standard 
non-standard 
non-standard 
non-standard 
non-standard 


Function Category Precision 


trig 

trig 

trig 

divide 
divide 
integer 
integer 
error function 
error function 
error function 
error function 
exp 

exp 

exp 

exp 

exp 

exp 

abs 

abs 

other 

other 
rounding 
rounding 
other 

other 
minmax 
minmax 
minmax 
minmax 
modrem 
modrem 
comparison 
comparison 
other 

other 

other 

other 

log 

log 

other 
rounding 
comparison 
comparison 
comparison 
comparison 


single 
double 
single 
double 
single 
word 
word 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
single 
single 
single 
double 
double 
single 
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Both 
SPU 
Both 
SPU 
Both 
Both 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
Both 
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Function Name 


isfinited2 
isfinitef4 
isgreaterd2 
isgreaterf4 
isgreaterequald2 
isgreaterequalf4 
islessd2 
islessf4 
islessequald2 
islessequalf4 
islessgreaterd2 
islessgreaterf4 
isunorderedd2 
isunorderedf4 
isinfd2 

isinff4 

isnand2 
isnanf4 
isnormald2 
isnormalf4 
Idexpd2 
Idexpf4 
Igammad2 
Igammaf4 
llabsi2 

IIdivi2 

lidivu2 

Ilrintd2 

Ilrintf4 
llroundd2 
llroundf4 
log10d2 

log 10f4 
log1pd2 
log1pf4 
log2d2 

log2f4 

logbd2 

logbf4 

logd2 

logf4 

modfd2 
modff4 
nearbyintd2 
nearbyintf4 


C99 Name 
isfinite 
isfinite 
isgreater 
isgreater 
isgreaterequal 
isgreaterequal 
isless 

isless 
islessequal 
islessequal 
islessgreater 
islessgreater 
isunordered 
isunordered 
isinf 

isinf 

isnan 

isnan 
isnormal 
isnormal 
Idexp 

Idexpf 
Igamma 
Igammaf 
llabs 

IIdiv 
non-standard 
Ilrint 

llrintf 

llround 
llroundf 
log10 

log10f 

log1p 

log1pf 

log2 

log2f 

logb 

logbf 

log 

logf 

modf 

modff 
nearbyint 
nearbyintf 


Function Category Precision 


comparison 
comparison 
comparison 
comparison 
comparison 
comparison 
comparison 
comparison 
comparison 
comparison 
comparison 
comparison 
comparison 
comparison 
comparison 
comparison 
comparison 
comparison 
comparison 
comparison 
other 

other 
gamma 
gamma 
integer 
integer 
integer 
rounding 
rounding 
rounding 
rounding 
log 

log 

log 

log 

log 


log 

modrem 
modrem 
rounding 
rounding 


double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double word 
double word 
double word 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 


SPU/PPU 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
SPU 
SPU 
SPU 
SPU 
SPU 
SPU 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
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Function Name 


negated2 
negatef4 
negatei4 
negatell2 
nextafterd2 
nextafterf4 
powd2 
powf4 
recipd2 
recipf4 
remainderd2 
remainderf4 
remquod2 
remquof4 
rintd2 
rintf4 
roundd2 
roundf4 
rsqrtd2 
rsqrtf4 
scalblind2 
scalbnf4 
signbitf4 
signbitd2 
sincosd2 
sincosf4 
sind2 

sinf4 
sinhd2 
sinhf4 
sqrtd2 
sqrtf4 
tand2 
tanf4 
tanhd2 
tanhf4 
tgammad2 
tgammaf4 
truncd2 
truncf4 


C99 Name 


non-standard 
non-standard 
non-standard 
non-standard 
nextafter 
nextafterf 
pow 

powf 
non-standard 
non-standard 
remainder 
remainderf 
remquo 
remquof 

rint 

rintf 

round 

roundf 
non-standard 
non-standard 
non-standard 
scalbnf 
signbit 
signbit 
non-standard 
non-standard 
sin 

sinf 

sinh 

sinhf 

sqrt 

sqrtf 

tan 

tanf 

tanh 

tanhf 
tgamma 
tgammaf 
trunc 

truncf 
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other 
other 
integer 
integer 
rounding 
rounding 
power 
power 
recip 
recip 
modrem 
modrem 
modrem 
modrem 
rounding 
rounding 
rounding 
rounding 
sqrt 

sqrt 
other 
other 
comparison 
comparison 
trig 

trig 

trig 

trig 

trig 

trig 

sqrt 

sqrt 

trig 

trig 

trig 

trig 
gamma 
gamma 
rounding 
rounding 


double 
single 
word 
double word 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
single 
double 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 
double 
single 


SPU/PPU 
SPU 
Both 
Both 
SPU 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
Both 
SPU 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
SPU 
Both 
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C. Special Cases 


Unless otherwise specified, each element of a SIMD result will adhere to either the C99 standard or the /EC 
60559:1989 standard. 


1. Rounding 


On the SPU, the full range of IEEE rounding modes is supported for double precision, but only round-toward-zero is 
supported for single precision. On the PPU, SIMD operations always use IEEE round-to-nearest mode. 


The mathematical accuracy of the SIMD functions assumes the default rounding mode. Accuracy may be 
compromised if the functions are called in another rounding mode. 


2. Special Operands 


On the PPU , NaN and Inf are recognized as special operands. 


On the SPU, all values passed to single-precision functions are treated as ordinary operands. NaN and Inf are not 
recognized as special single-precision operands; however, they are recognized as special double-precision operands, 
and SIMD functions check for them, as described in C99, IEC 60559:1989 and the SIMD function specifications. See 
Synergistic Processor Unit Instruction Set Architecture for details. 


On both the PPU and the SPU, single precision floating-point denormal inputs are coerced to zero unless otherwise 
noted. 


3. Error Conditions 


A domain error occurs if an input argument is outside the domain over which the mathematical function is defined. 
The description of each function lists any required domain errors. The resulting vector element is undefined for all 
corresponding element input arguments which contain a domain error and no exception or error is reported. 


A range error occurs when the mathematical result cannot be represented in an object of the specified type. When a 
range error occurs, the resulting vector element is either HUGE VAL (for double precision results) or HUGE VALF (for 
single precision results). Integer arithmetic function results are undefined when they cannot be represented. 








4. Exceptions 


The SIMD library functions have an undefined effect on the exception flags in the SPU floating-point status and 
control register (FPSCR). SPU functions on double precision arguments set exception bits in the FPSCR that can be 
tested by calling the routine fegetexcept, as documented in fenv.h. 


The SPU does not raise hardware traps for single-precision fp exceptions; PPU SIMD operations do have hardware 
support for a subset of the C99 fp exceptions. 
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2. SIMD Function Specifications 


This chapter contains descriptions of the SIMD math functions, their arguments, and their return values. Where 
necessary, accuracy information is provided to clarify expected specific behavior. All functions are available on both 
the PPU and SPU unless otherwise noted. 


A. Type Definitions 


The following type definitions are used for function return values: 


divi4_t: Remainder/Quotient Struct for Vector Signed Int 


typedef struct divi4 s ( 
vector signed int quot; 
vector signed int rem; 
) divi4 t; 


Structures of this type are used to hold the return value of divi4(). The member quot contains the quotient and the 
member rem contains the remainder of the division. 


divu4 t: Remainder/Quotient Struct for Vector Unsigned Int 
typedef struct divu4 s { 
vector unsigned int quot; 
vector unsigned int rem; 
) divu4 t; 


Structures of this type are used to hold the return value of divu4 ().The member quot contains the quotient and the 
member rem contains the remainder of the division. 


IIdivi2_t: Remainder/Quotient Struct for Vector Signed Long Long (SPU Only) 
typedef struct lldivi2 s ( 


vector signed long long quot; 

vector signed long long rem; 
) lldivi2 t; 
Structures of this type are used to hold the return value of 11divi2 () .The member quot contains the quotient and 
the member rem contains the remainder of the division. 


lldivu2 t: Remainder/Quotient Struct for Vector Unsigned Long Long (SPU Only) 


typedef struct lldivu2 s ( 
vector unsigned long long quot; 
vector unsigned long long rem; 
} ll1divu2 t; 


Structures of this type are used to hold the return value of 11divu2 () .The member quot contains the quotient and 


the member rem contains the remainder of the division. 


Ilroundf4_t: Vector of Four Long Long (SPU Only) 


typedef struct llroundf4 s ( 
vector signed long long vll[2]; 
) llroundf4 t; 


Structures of this type are used to hold signed long long data corresponding to a vector of four elements. 
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B. Function Descriptions 
In the function descriptions that follow, a subscript is used to indicate a vector element. For example, element / of 


vector x is shown as x;. 


absi4: Absolute Value of Integer 


(vector signed int) absi4 (vector signed int x) 
A vector signed int is returned that contains the absolute value of each corresponding element of vector signed int x. 


If the absolute value of x; cannot be represented, the corresponding result is undefined and no error is reported. 


acosd2: Arccosine of Double (SPU Only) 


(vector double) acosd2 (vector double x); 


A vector double is returned that contains the angles whose cosines correspond to the respective elements in vector 
double x. Each element in the result is in the range [0, pi] radians. 


If the absolute value of x; is greater than 1, the corresponding result is undefined and no error is reported. 


acosf4: Arccosine of Float 


(vector float) acosf4 (vector float x); 


A vector float is returned that contains the angles whose cosines correspond to the respective elements in vector float 
x. Each element in the result is in the range [0, pi] radians. 


If the absolute value of x; is greater than 1, the corresponding result is undefined and no error is reported. 


acoshd2: Hyperbolic Arccosine of Double (SPU Only) 


(vector double) acoshd2 (vector double x) 


A vector double is returned that contains the nonnegative hyperbolic arccosines of the corresponding elements of 
vector double x. 


If the value of x; is less than 1, the corresponding result is undefined and no error is reported. 


acoshf4: Hyperbolic Arccosine of Float 


(vector float) acoshf4 (vector float x) 


A vector float is returned that contains the nonnegative hyperbolic arccosines of the corresponding elements of vector 
float x. 


If the value of x; is less than 1, the corresponding result is undefined and no error is reported. 


asind2: Arcsine of Double (SPU Only) 


(vector double) asind2 (vector double x); 


A vector double is returned that contains the angles whose sines correspond to the respective elements in vector 
double x. Each element in the result is in the range [-pi/2, *pi/2] radians. 


If the absolute value of x; is greater than 1, the corresponding result is undefined and no error is reported. 


asinf4: Arcsine of Float 


(vector float) asinf4 (vector float x); 


A vector float is returned that contains the angles whose sines correspond to the respective elements in vector float x. 
Each element in the result is in the range [-pi/2, *pi/2] radians. 


If the absolute value of x, is greater than 1, the corresponding result is undefined and no error is reported. 
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asinhd2: Hyperbolic Arcsine of Double (SPU Only) 
(vector double) asinhd2 (vector double x) 
A vector double is returned that contains the nonnegative hyperbolic arcsines of the corresponding elements of vector 
double x. 


asinhf4: Hyperbolic Arcsine of Float 
(vector float) asinhf4 (vector float x) 
A vector float is returned that contains the nonnegative hyperbolic arcsines of the corresponding elements of vector 
float x. 


atand2: Tangent of Double (SPU Only) 
(vector double) atand2 (vector double x); 
A vector double is returned that contains the angles whose tangents correspond to the respective elements of vector 


double x. Each element in the result is in the range [-pi/2, *pi/2] radians. 


atanf4: Tangent of Float 
(vector float) atanf4 (vector float x); 
A vector float is returned that contains the angles whose tangents correspond to the respective elements of vector 


float x. Each element in the result is in the range [-pi/2, *pi/2] radians. 


atanhd2: Hyperbolic Arctangent of Double (SPU Only) 


(vector double) atanhd2 (vector double x) 


A vector double is returned that contains the nonnegative hyperbolic arctangents of the corresponding elements of 
vector double x. 


If the absolute value of x; is greater than 1, the corresponding result is undefined and no error is reported. 


atanhf4: Hyperbolic Arctangent of Float 


(vector float) atanhf4 (vector float x) 


A vector float is returned that contains the nonnegative hyperbolic arctangents of the corresponding elements of 
vector float x. 


If the absolute value if x, is greater than 1, the corresponding result is undefined and no error is reported. 


On the SPU, if the absolute value of x, is equal to 1, the corresponding element of the result will be returned as 
HUGE VALF and no erroris reported. 





atan2d2: Arctangent of Double Quotient (SPU Only) 


(vector double) atan2d2 (vector double y, vector double x); 


A vector double is returned that contains the angles whose tangents are y;/ x; for corresponding elements of vector 
double y and vector double x. Each element in the result is within the range [-pi, +pi] radians. 


If x; and y; are zero, the corresponding element of the result is undefined and no error is reported. 


atan2f4: Arctangent of Float Quotient 


(vector float) atan2f4 (vector float y, vector float x); 


A vector float is returned that contains the angles whose tangents are y;/ x; for corresponding elements of vector 
float y and vector float x. Each element in the result is within the range [-pi, *pi] radians. 


If x; and y, are zero, the corresponding element of the result is undefined, and no error is reported. 
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cbrtd2: Cube Root of Double (SPU Only) 
(vector double) cbrtd2 (vector double x); 


143 


A vector double is returned that contains the real cube roots, x; ^" of the corresponding elements of vector double x. 


cbrtf4: Cube Root of Float 
(vector float) cbrtf4 (vector float x); 


1/3 


A vector float is returned that contains the real cube roots, x; ^" of the corresponding elements of vector float x. 


ceild2: Ceiling of Double (SPU Only) 
(vector double) ceild2 (vector double x); 
A vector double is returned that contains the smallest integer values, expressed as floating-point numbers, that are 


not less than the corresponding elements of vector double x. 


ceilf4: Ceiling of Float 
(vector float) ceilf4 (vector float x); 
A vector float is returned that contains the smallest integer values, expressed as floating-point numbers, that are not 


less than the corresponding elements of vector float x. 


copysignd2: Copy Sign of Double (SPU Only) 
(vector double) copysignd2 (vector double x, vector double y); 
A vector double is returned that contains the magnitude of the corresponding element of vector double x and the sign 


of the corresponding element of vector double y. 


copysignf4: Copy Sign of Float 
(vector float) copysignf4 (vector float x, vector float y); 
A vector float is returned that contains the magnitude of the corresponding element of vector float x and the sign of 


the corresponding element of vector float y. 


cosd2: Cosine of Double (SPU Only) 
(vector double) cosd2 (vector double x); 
A vector double is returned that contains the cosines of the corresponding elements of vector double x. 
The results of cosd2 () may not be accurate for very large values of x, but no error is reported. Implementations 


should document the point at which accuracy is lost. 


cosf4: Cosine of Float 
(vector float) cosf4 (vector float x); 
A vector float is returned that contains the cosines of the corresponding elements of vector float x. 
The results of cos£4 () are may not be accurate for very large values, but no error is reported. Implementations 


should document the point at which accuracy is lost. 


coshd2: Hyperbolic Cosine of Double (SPU Only) 


(vector double) coshd2 (vector double x) 


A vector double is returned that contains the hyperbolic cosines of the corresponding elements of vector double x. 
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coshf4: Hyperbolic Cosine of Float 


(vector float) coshf4 (vector float x) 
A vector float is returned that contains the hyperbolic cosines of the corresponding elements of vector float x. 


On the SPU, element values of the result that are greater than HUGE VALF are returned as HUGE VALF, and no error 
is reported. 








divd2: Divide Doubles (SPU Only) 
(vector double) divd2 (vector double x, vector double y); 


A vector double is returned that contains the quotient x; / y;, for the corresponding elements of vector double x and 
vector double y. This function handles special cases as follows: 


e |f either input is NaN, the result is NaN 

e For Inf/Inf or 0/0, the result is NaN 

e For finite/O, the result is Inf with sign = sign(x)/sign( y) 
e For finite/t Inf, the result is O with sign = sign(x)/sign( y) 


divf4: Divide Floats 
(vector float) divf4 (vector float x, vector float y); 
A vector float is returned that contains the quotients x; /y;, for the corresponding elements of vector float x and 
vector float y. This function handles special cases as follows: 
e |f either input is NaN, the result is NaN 
e For Inf/Inf or 0/0, the result is NaN 
e For finite/O, the result is Inf with sign = sign(x)/sign( y) 
e  Forfinite/£Inf, the result is O with sign = sign(x)/sign( y) 





On the SPU, if y, is zero, the result is HUGE VALF with sign = sign(x)/sign( y). 


divi4: Divide Integer 
(divi4 t) divi4 (vector signed int x, vector signed int y) 


Each element of vector signed int x is divided by the corresponding element of vector signed int y, and the result is 
returned in a structure of type divi4 t() that contains a vector of corresponding quotients and a vector of 
corresponding remainders. 


Each element in the structure member quot is the algebraic quotient truncated towards zero. Each element in the 
structure member remis the corresponding remainder, such that x; == quot * y; + rem. 


If y; is zero, the corresponding element of the resulting quotient is zero. 


divu4: Divide Unsigned Integer 


(divu4 t) divu4 (vector unsigned int x, vector unsigned int y) 


Each element of vector unsigned int x is divided by the corresponding element of vector unsigned int y, and the result 
is returned in a structure of type divu4 t() that contains a vector of corresponding quotients and a vector of 
corresponding remainders. 


Each element in the structure member quot is the algebraic quotient truncated towards zero. Each element in the 
structure member remis the corresponding remainder, such that x; == quot * y; + rem. 


If y; is zero, the corresponding element of the resulting quotient is zero. 
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erfcd2: Complementary Error Function Double (SPU Only) 
(vector double) erfcd2 (vector double x) 
A vector double is returned that contains the complementary error functions of the corresponding elements of vector 
double x. 


erfcf4: Complementary Error Function Float 
(vector float) erfcf4 (vector float x) 
A vector float is returned that contains the complementary error functions of the corresponding elements of vector 
float x. 


erfd2: Error Function Double (SPU Only) 


(vector double) erfd2 (vector double x) 


A vector double is returned that contains the error functions of the corresponding elements of vector double x. 


erff4: Error Function Float 


(vector float) erff4 (vector float x) 


A vector float is returned that contains the error functions of the corresponding elements of vector float x. 


expd2: e Raised to the Power of Double (SPU Only) 


(vector double) expd2 (vector double x); 


A vector double is returned that contains the corresponding exponentials e*, for each element of vector double x. 


expf4: e Raised to the Power of Float 
(vector float) expf4 (vector float x); 
A vector float is returned that contains the corresponding exponentials e*, for each element of vector float x. 


On the SPU, element values of the result that are greater than HUGE VALF are returned as HUGE VALF, and no error 
is reported. 








exp2d2: 2 Raised to the Power of Double (SPU Only) 


(vector double) exp2d2 (vector double x); 


A vector double is returned that contains the corresponding exponentials 2*,, for each element of vector double x. 


exp2f4: 2 Raised to the Power of Float 
(vector float) exp2f4 (vector float x); 
A vector float is returned that contains the corresponding exponentials 2*;, for each element of vector float x. 


On the SPU, element values of the result that are greater than HUGE VALF are returned as HUGE VALF and no error 
is reported. 








expm1d2: e Raised to the Power of Double Minus 1 (SPU Only) 


(vector double) expmld2 (vector double x); 


A vector double is returned that contains the exponential minus 1, e*; — 1, for corresponding elements of vector 
double x. 


This function returns mathematically accurate values, even when x; is near 0, or when exp (x;)-1.0 would return 
bad values due to floating-point cancellation errors. 
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expm1f4: e Raised to the Power of Float Minus 1 
(vector float) expmlf4 (vector float x); 
A vector float is returned that contains the exponential minus 1, e*; — 1, for corresponding elements of vector float x. 
This function returns mathematically accurate values, even when an element of x, is near 0, or when expf (xi) - 


1. 0f would return bad values due to floating-point cancellation errors. 


fabsd2: Absolute Value Double (SPU Only) 
(vector double) fabsd2 (vector double x); 


A vector double is returned that contains the absolute values, | x, |, for corresponding elements of vector double x. 


fabsf4: Absolute Value Float 
(vector float) fabsf4 (vector float x); 


A vector float is returned that contains the absolute values, | x; |, for corresponding elements of vector float x. 


fdimd2: Subtract Staying Non-Negative Double (SPU Only) 
(vector double) fdimd2 (vector double x, vector double y); 
A vector double is returned that contains the larger of (x;-y;) and zero, for corresponding elements of vector double 


x and vector double y. 


fdimf4: Subtract Staying Non-Negative Float 
(vector float) fdimf4 (vector float x, vector float y); 
A vector float is returned that contains the larger of (x;-y;) and zero, for corresponding elements of vector float x 


and vector float y. 


floord2: Floor Double (SPU Only) 
(vector double) floord2 (vector double x); 
A vector double is returned that contains the largest integer values, expressed as floating-point numbers, that are not 


greater than the corresponding elements of vector double x. 


floorf4: Floor Float 
(vector float) floorf4 (vector float x); 
A vector float is returned that contains the largest integer values, expressed as floating-point numbers, that are not 


greater than the corresponding elements of vector float x. 


fmad2: Fused Multiply and Add Double (SPU Only) 
(vector double) fmad2 (vector double x, vector double y, vector double z); 
A vector double is returned that contains the results of the calculation of (x; * y; + z;), for the corresponding 


elements of vector double x, vector double y, and vector double z. Intermediate results are of arbitrary precision. 


fmaf4: Fused Multiply and Add Float 


(vector float) fmaf4 (vector float x, vector float y, vector float z); 


A vector float is returned that contains the results of the calculation of (x; * y; + z;), for the corresponding elements 
of vector float x, vector float y, and vector float z. Intermediate results are of arbitrary precision. 
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fmaxd2: Maximum Double (SPU Only) 


(vector double) fmaxd2 (vector double x, vector double y); 
A vector double is returned that contains the larger (more positive) of x, and y,, for corresponding elements of vector 
double x and vector double y. 


fmaxf4: Maximum Float 


(vector float) fmaxf4 (vector float x, vector float y); 


A vector float is returned that contains the larger (more positive) of x; and y;, for corresponding elements of vector 
float x and vector float y. 


On the SPU, this function does not coerce denormals to zero. Instead, it compares them as normal values even 
though the SPU's floating-point instructions do not. 
fmind2: Minimum Double (SPU Only) 


(vector double) fmind2 (vector double x, vector double y); 
A vector double is returned that contains the smaller (more negative) of x; and y,, for corresponding elements of 
vector double x and vector double y. 


fminf4: Minimum Float 


(vector float) fminf4 (vector float x, vector float y); 


A vector float is returned that contains the smaller (more negative) of x; and y;, for corresponding elements of vector 
float x and vector float y. 


On the SPU, this function does not coerce denormals to zero. Instead, it compares them as normal values even 
though the SPU's floating-point instructions do not. 


fmodd2: Modulus Double (SPU Only) 


(vector double) fmodd2 (vector double x, vector double y); 


A vector double is returned that contains the remainder of x;/ y;, for corresponding elements of vector double x and 
vector double y, as described below: 


e |fy;is 0, the result is 0 
e Otherwise, the function determines the unique signed integer value i such that the returned element is 
xi- i * y, With the same sign as x; and magnitude less than | y, | 
fmodf4: Modulus Float 


(vector float) fmodf4 (vector float x, vector float y); 


A vector float is returned that contains the remainder of x;/ y;, for corresponding elements of vector float x and 
vector float y, as defined below: 


e If y, is 0, the result is 0 
e Otherwise, fmodf4 () determines the unique signed integer value i such that the returned element is 
x;-i* y, with the same sign as x; and magnitude less than | y, | 
fpclassifyd2: Classify Double (SPU Only) 
(vector signed long long) fpclassifyd2 (vector double x) 
A vector signed long long is returned that contains the floating-point classifications for corresponding elements of 


vector double x. The classifications, which are defined in math.h, are FP NAN, FP. INFINITE, FP NORMAL, 
FP SUBNORMAL, and FP ZERO. 
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fpclassifyf4: Classify Float 
(vector signed int) fpclassifyf4 (vector float x) 
A vector signed int is returned that contains the floating-point classifications for corresponding elements of vector float 


x. The classifications, which are defined in math.h, are FP_NAN, FP INFINITE, FP NORMAL, FP_ SUBNORMAL, and 
FP ZERO. 











On the SPU, the resulting vector will never contain FP_NAN or FP_INFINITE. 


frexpd2: Represent Double as Fraction and Exponent (SPU Only) 
(vector double) frexpd2 (vector double x, vector signed long long *pexp); 
A vector double is returned that contains normalized fractions, and a vector signed long long is stored in *pexp that 


contains exponent integers. Each fraction element £rac and each exponent integer element exp represents the 
value of the corresponding element of x, such that: 


e Every element of | £rac| is in the interval [1/2, 1) or is zero 
e x;-- frac * 277? 

e Ifx;is 0, the corresponding element of *pexp is also zero 

e If x; is NaN, the corresponding result is NaN and the corresponding element of *pexp is undefined. 


e = If x; is infinite, the corresponding result is infinite and the corresponding element of *pexp is 
undefined. 
frexpf4: Represent Float as Fraction and Exponent 


(vector float) frexpf4 (vector float x, vector signed int *pexp); 


A vector float is returned that contains normalized fractions, and a vector signed int is stored in *pexp that contains 
exponent integers. Each fraction element frac and each exponent integer element exp represents the value of the 
corresponding element of x, such that: 


e Every element of | frac| is in the interval [1/2, 1) or is zero 

e x;-- frac * 277? 

e Ifx;is 0, the corresponding element of *pexp is also zero 

e If x; is NaN, the corresponding result is NaN and the corresponding element of *pexp is undefined. 


e = If x; is infinite, the corresponding result is infinite and the corresponding element of *pexp is 
undefined. 


hypotd2: Hypotenuse Double (SPU Only) 
(vector double) hypotd2 (vector double x, vector double y); 
A vector double is returned that contains the square root of x^; + y^; without undue overflow or underflow, for 


corresponding elements of vector double x and vector double y. 


hypotf4: Hypotenuse Float 
(vector float) hypotf4 (vector float x, vector float y); 
A vector float is returned that contains the square root of x°; + y^; without undue overflow or underflow, for 


corresponding elements of vector float x and vector float y. 


ilogbd2: Integer Exponent of Double (SPU Only) 
(vector signed long long) ilogbd2 (vector double x); 


A vector signed long long is returned that contains the elements defined below, for corresponding elements of vector 
double x. 


e If x; is not-a-number (NaN), the value is the macro FP. ILOGBNAN 
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e If x; is equal to zero, the value is the macro FP ILOGBO 
e If x, is equal to positive or negative Inf, the value is the macro FP. ILOGBO 
e Otherwise, the result is (int)logb(xi) 
ilogbf4: Integer Exponent of Float 
(vector signed int) ilogbf4 (vector float x); 
A vector signed int is returned that contains the elements defined below, for corresponding elements of vector float x. 


e If x; is not-a-number (NaN), the value is the macro FP. ILOGBNAN 

e If x; is equal to 0, the value is the macro FP ILOGBO 

e If x; is equal to positive or negative Inf, the value is the macro FP. ILOGBO 
e Otherwise, the result is (int) logb (xi) 


Because the SPU treats single-precision Inf and NaN codes as regular floating-point numbers, ilogbf4 returns a 
result of 128 for these numbers. For compatibility with the double function ilogb(), FP ILOGBNAN is set to 
INT MAX. 


irintf4: Nearest Integer Float 


(vector signed int) irintf4 (vector float x) 


A vector signed int is returned that contains the nearest integer to the corresponding element of vector float x, 


consistent with the current rounding mode. If the rounded value is outside the range of the return type, the numeric 
result is unspecified. 


On the SPU, the rounding mode for floats is always towards zero. 


iroundf4: Round Float to Nearest Integer 


(vector signed int) iroundf4 (vector float x) 
A vector signed int is returned that contains the rounded integer value of the corresponding element of vector float x. 


Elements are rounded to the nearest value; halfway cases are rounded away from zero, regardless of the current 
rounding direction. 


If the rounded value is outside the range of the return type, the numeric result is unspecified. 


is0denormd2: 0 or Denormalized Double (SPU Only) 


(vector unsigned long long) isOdenormd2 (vector double x); 


A vector unsigned long long is returned that contains the elements defined below for corresponding elements of 
vector double x. 


e All bits of the resulting element are set to 1 if x; is a denormalized value or zero 
e Zero otherwise 
isOdenormf4: 0 or Denormalized Float 
(vector unsigned int) isOdenormf4 (vector float x); 


A vector unsigned int is returned that contains the elements defined below for corresponding elements of 
vector float x. 


e All bits of the resulting element are set to 1 if x; is a denormalized value or zero 
e Zero otherwise 
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isequald2: Compare Equal Double (SPU Only) 
(vector unsigned long long) isequald2 (vector double x, vector double y); 


A vector unsigned long long is returned that contains the elements defined below for corresponding elements of 
vector double x and vector double y. 


e  Allbits of the resulting element are set to 1 if x, and y, are equal 
e Zero otherwise 


The function correctly compares denormalized numbers. If either input is NaN, the comparison result is false (zero). If 
both inputs are Inf with same sign, the inputs are considered equal. The values 0 and -0 are considered equal. 


isequalf4: Compare Equal Float 


(vector unsigned int) isequalf4 (vector float x, vector float y); 


A vector unsigned int is returned that contains the elements defined below for corresponding elements of vector float 
x and vectory float y. 


e All bits of the resulting element are set to 1 if x, and y, are equal 
e Zero otherwise 


The function correctly compares denormalized numbers. If either input is NaN, the comparison is false (zero). If both 
inputs are Inf with same sign, the inputs are considered equal. The values 0 and -0 are considered equal. 


isfinited2: Double is Finite (SPU Only) 


(vector unsigned long long) isfinited2 (vector double x) 


A vector unsigned long long is returned that contains the elements defined below for corresponding elements of 
vector double x. 


e All bits of the resulting element are set to 1 if x, is finite 


e Zero otherwise 


isfinitef4: Float is Finite 


(vector unsigned int) isfinitef4 (vector float x) 
A vector unsigned int is returned that contains elements defined below for corresponding elements of vector float x. 


e All bits of the resulting element are set to 1 if x, is finite 
e Zero otherwise 


On the SPU, infinite values are not representable in single precision. Therefore, all bits of the resulting element are 
set to 1 regardless of the value of x;. 


isgreaterequald2: Greater or Equal Double (SPU Only) 
(vector unsigned long long) isgreaterequald2 (vector double x, vector double y); 


A vector unsigned long long is returned that contains elements defined below for corresponding elements of vector 
double x and vector double y. 


e All bits of the resulting element are set to 1 if x; equals to or greater than y; 
e Zero otherwise 


The function correctly compares denormalized numbers. If either element of the input is NaN, the comparison is false. 
If both elements of the inputs are Inf with the same sign, the inputs are considered equal. The values 0 and -0 are 
considered equal. 
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isgreaterequalf4: Greater or Equal Float 


(vector unsigned int) isgreaterequalf4 (vector float x, vector float y); 


A vector unsigned int is returned that contains elements defined below for corresponding elements of vector float x 
and vector float y. 


e  Allbits of the resulting element are set to 1 if x; equals to or greater than y; 
e Zero otherwise 


The function correctly compares denormalized numbers. If either element of the input is NaN, the comparison is false. 
If both elements of the inputs are Inf with the same sign, the inputs are considered equal. The values 0 and -0 are 
considered equal. 


isgreaterd2: Greater than Double (SPU Only) 


(vector unsigned long long) isgreaterd2 (vector double x, vector double y); 


A vector unsigned long long is returned that contains elements defined below for corresponding elements of vector 
double x and vector double y. 


e All bits of the resulting element are set to 1 if x, greater than y; 


e Zero otherwise 


The function correctly compares denormalized numbers. 


isgreaterf4: Greater than Float 


(vector unsigned int) isgreaterf4 (vector float x, vector float y); 


A vector unsigned int is returned that contains elements defined below for corresponding elements of vector float x 
and vector float y. 


e All bits of the resulting element are set to 1 if x; greater than y; 
e Zero otherwise 


The function correctly compares denormalized numbers. 


isinfd2: Double is Infinity (SPU Only) 


(vector unsigned long long) isinfd2 (vector double x) 


A vector unsigned long long is returned that contains elements defined below, for corresponding elements of vector 
double x. 


e All bits of the resulting element are set to 1 if x; is infinite 


e Zero otherwise 


isinff4: Float is Infinity 
(vector unsigned int) isinff4 (vector float x) 


A vector unsigned long int is returned that contains elements defined below for corresponding elements of vector 
float x. 


e All bits of the resulting element are set to 1 if x, is infinite 
e Zero otherwise 


On the SPU, infinite values are not representable in single precision. Therefore, all bits of the resulting element are 
set to zero, regardless of the value of x;. 
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islessd2: Double is Less Than (SPU Only) 


(vector unsigned long long) islessd2 (vector double x, vector double y) 


A vector unsigned long long is returned that contains elements defined below, for corresponding elements of vector 
double x and vector double y. 


e All bits of the resulting element are set to 1 if x, less than y; 


e Zero otherwise 


The function correctly compares denormalized numbers. 


islessequald2: Double is Less Than or Equal To (SPU Only) 


(vector unsigned long long) islessequald2 (vector double x, vector double y) 


A vector unsigned long long is returned that contains elements defined below for corresponding elements of vector 
double x and vector double y. 


e All bits of the resulting element are set to 1 if x, less than or equal y; 
e Zero otherwise 


The function correctly compares denormalized numbers. 


islessequalf4: Float is Less Than or Equal To 


(vector unsigned int) islessequalf4 (vector float x, vector float y) 


A vector unsigned int is returned that contains elements defined below for corresponding elements of vector float x 
and vector float y. 


e All bits of the resulting element are set to 1 if x; less than or equal y; 
e Zero otherwise 


The function correctly compares denormalized numbers. 


islessf4: Float is Less Than 


(vector unsigned int) islessf4 (vector float x, vector float y) 


A vector unsigned int is returned that contains elements defined below for corresponding elements of vector float x 
and vector float y. 


e All bits of the resulting element are set to 1 if x; less than y; 
e Zero otherwise 


The function correctly compares denormalized numbers. 


islessgreaterd2: Double is Less Than or Greater Than (SPU Only) 





(vector unsigned long long) islessgreaterd2 (vector double x, vector double y) 


A vector unsigned long long is returned that contains elements defined below for corresponding elements of vector 
double x and vector double y. 


e All bits of the resulting element are set to 1 if x; less than or greater than y; 
e Zero otherwise 


The function correctly compares denormalized numbers. 
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islessgreaterf4: Float is Less Than or Greater Than 


(vector unsigned int) islessgreaterf4 (vector float x, vector float y) 





A vector unsigned int is returned that contains elements defined below, for corresponding elements of vector float x 
and vector float y. 


e All bits of the resulting element are set to 1 if x; less than or greater than y; 
e Zero otherwise 


The function correctly compares denormalized numbers. 


isnand2: Double is NaN (SPU Only) 


(vector unsigned long long) isnand2 (vector double x) 


A vector unsigned long long is returned that contains elements defined below for corresponding elements of vector 
double x. 


e All bits of the resulting element are set to 1 if x; is a NaN 


e Zero otherwise 


isnanf4: Float is NaN 
(vector unsigned int) isnanf4 (vector float x) 
A vector unsigned int is returned that contains elements defined below, for corresponding elements of vector float x. 
e All bits of the resulting element are set to 1 if x; is a NaN 


e Zero otherwise 


On the SPU, NaN is not representable in single precision. Therefore, all bits of the resulting element are set to zero, 
regardless of the value of x;. 


isnormald2: Double is Normal (SPU Only) 


(vector unsigned long long) isnormald2 (vector double x) 


A vector unsigned long long is returned that contains elements defined below for corresponding elements of vector 
double x. 


e All bits of the resulting element are set to 1 if x; is normal, not a Nan, or an infinity 
e Zero otherwise 
isnormalf4: Float is Normal 
(vector unsigned int) isnormalf4 (vector float x) 
A vector unsigned int is returned that contains elements defined below for corresponding elements of vector float x. 
e all bits of the resulting element are set to 1 if x; is normal, not a NaN, or an infinity 
e zero otherwise 
isunorderedd2: Double is Unordered (SPU Only) 
(vector unsigned long long) isunorderedd2 (vector double x, vector double y) 


A vector unsigned long long is returned that contains elements defined below for corresponding elements of vector 
double x and vector double y. 


e All bits of the resulting element are set to 1 if x; is unordered to the element of y, 
e Zero otherwise 
NaN is unordered to any operand, including NaN itself. 
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isunorderedf4: Float is Unordered 


(vector unsigned int) isunorderedf4 (vector float x, vector float y) 


A vector unsigned int is returned that contains elements defined below for corresponding elements of vector float x 
and vector float y. 

e  Allbits of the resulting element are set to 1 if x; is unordered to y; 

e Zero otherwise 


NaN is unordered to any operand, including NaN itself. On the SPU, NaN does not exist in single precision. Therefore, 
this function will always return zero. 


Idexpd2: Multiply Double by 2 Raised to its Power (SPU Only) 


(vector double) ldexpd2 (vector double x, vector signed long long ex); 


A vector double is returned that contains x, * 2**, for the corresponding elements of vector double x and vector 
signed long long ex. For large elements of ex; (overflow), the element in the result saturates to HUGE VAL with an 
appropriate sign. For small elements of ex (underflow), the corresponding element of the result is 0. 





Idexpf4: Multiply Float by 2 Raised to its Power 
(vector float) ldexpf4 (vector float x, vector signed int ex); 
A vector float is returned that contains x, * 2**, for the corresponding elements of vector float x and vector signed 


int ex. For large elements of ex (overflow), the element in the result saturates to HUGE VALF with an appropriate sign. 
For small ex; (underflow), the corresponding element of the result is 0. 





Igammad2: Natural Log of Gamma Function of Double (SPU Only) 
(vector double) lgammad2 (vector double x) 
A vector double is returned that contains the natural logarithm of the absolute value of the result of the gamma 


function for the corresponding elements of vector double x. 


Igammaf4: Natural Log of Gamma Function of Float 
(vector float) lgammaf4 (vector float x) 
A vector float is returned that contains the natural logarithm of the absolute value of the result of the gamma function 


for the corresponding element of vector float x. 


llabsi2: Absolute Value Long Long (SPU Only) 


(vector long long) llabsi2 (vector signed long long x) 


A vector long long is returned that contains the absolute value, | x;| of the corresponding element of vector signed 
long long x. 


If the absolute value of x, cannot be represented, the corresponding result is undefined and no error is reported. 


IIdivi2: Divide Long Long (SPU Only) 
(lldivi2 t) lldivi2 (vector signed long long x, vector signed long long y) 


Each element of vector signed long long x is divided by each element of vector signed long long y, and the result is 
returned in a structure of type 11divi2 t (), which contains a vector of quotients and a vector of remainders. 


Each element of the vector in the structure member quot is the algebraic quotient truncated towards zero. Each 
element of the vector in the structure member remis the corresponding remainder, such that x; == quot * y; + rem. 


If y; is zero, the corresponding element of the resulting quotient is zero. 
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lIldivu2: Divide Unsigned Long Long (SPU Only) 


(lldivu2 t) lldivu2 (vector unsigned long long x, vector unsigned long long y) 


Each element of vector unsigned long long x is divided by each element of vector unsigned long long y, and the 
result is returned in a structure of type 11divu2_t (), containing a vector of quotients and a vector of remainders. 


Each element of the vector in the structure member quot is the algebraic quotient truncated towards zero. Each 
element of the vector in the structure member rem is the corresponding remainder, such that x; == quot * y; + rem. 


If y, is zero, the corresponding element of the resulting quotient is zero. 


llrintd2: Find Nearest Long Long of Double (SPU Only) 


(vector signed long long) llrintd2 (vector double x) 


A vector signed long long is returned that contains the nearest long long integer to the corresponding element of 
vector double x consistent with the current rounding mode. If the rounded value is outside the range of the return type, 
the numeric result is unspecified. 


llrintf4: Find Nearest Long Long of Float (SPU Only) 
(llroundf4 t) llrintf4 (vector float x) 


A structure of type 11roundf4 t() is returned that contains the nearest long long integer to the corresponding 
element of vector float x consistent with the current rounding mode. If the rounded value is outside the range of the 
return type, the numeric result is unspecified. 


On the SPU the rounding mode for floats is always towards zero. 


Ilroundd2: Round Double to Nearest Long Long (SPU Only) 


(vector signed long long) llroundd2 (vector double x) 


A vector signed long long is returned that contains the corresponding elements of vector double x rounded to the 
nearest value, rounding halfway values away from 0 regardless of the current rounding direction. If the rounded value 
is outside the range of the return type, the numeric result is unspecified. 


lIroundf4: Round Float to Nearest Long Long (SPU Only) 
(llroundf4 t) llroundf4 (vector float x) 


A structure of type 11roundf4 t() is returned that contains the corresponding elements of vector float x rounded to 
the nearest value, rounding halfway cases away from 0 regardless of the current rounding direction. If the rounded 
value is outside the range of the return type, the numeric result is unspecified. 


logd2: Natural Log of Double (SPU Only) 
(vector double) logd2 (vector double x); 


A vector double is returned that contains the natural logarithms of the corresponding elements of vector double x. 


If x, is negative, the corresponding result is undefined and no error is reported. 


logf4: Natural Log of Float 
(vector float) logf4 (vector float x); 


A vector float is returned that contains the natural logarithms of the corresponding elements of vector float x. 
If x, is negative, the corresponding result is undefined and no error is reported. 


If x; is zero, the result is -HUGE_VALF. 
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log10d2: Log Base 10 of Double (SPU Only) 
(vector double) logl0d2 (vector double x); 
A vector double is returned that contains the base-10 logarithm of the corresponding elements of vector double x. 


If x; is negative, the corresponding result is undefined and no error is reported. 
log10f4: Log Base 10 of Float 
(vector float) log10f4 (vector float x); 


A vector float is returned that contains the base-10 logarithm of the corresponding elements of vector float x. 


If x, is negative, the corresponding result is undefined and no error is reported. 





If x; is zero, the result is -HUGE VALF. 


log1pd2: Natural Log of Double Plus 1 (SPU Only) 


(vector double) loglpd2 (vector double x); 


A vector double is returned that contains the natural logarithm of 1 + x; for the corresponding elements of vector 
double x. 


The function returns mathematically accurate values even when the corresponding element of x; is near zero. 


If x; is less than -1, the corresponding result is undefined and no error is reported. 


log1pf4: Natural Log of Float Plus 1 
(vector float) loglpf4 (vector float x); 
A vector float is returned that contains the natural logarithms of 1 + x;, for corresponding elements of vector float x. 


The function returns mathematically accurate values even when the corresponding element of x; is near zero. If an 
element of x; is -1, the resultis -HUGE VALF. 





If x; is less than -1, the corresponding result is undefined and no error is reported. 
log2d2: Log Base 2 of Double (SPU Only) 
(vector double) log2d2 (vector double x); 
A vector double is returned that contains the base-2 logarithm of the corresponding elements of vector double x. 
If x; is less than 0, the corresponding result is undefined and no error is reported. 
log2f4: Log Base 2 of Float 
(vector float) log2f4 (vector float x); 


A vector float is returned that contains the base-2 logarithm of the corresponding elements of vector float x. 


If x; is zero, the result is -HUGE VALF. 





If x; is less than zero, the corresponding result is undefined and no error is reported. 
logbd2: Represent Double as Fraction Greater Than 1 and Exponent (SPU Only) 
(vector double) logbd2 (vector double x); 


An integer exponent ex; and a fraction frac, that represents the value of a finite element are determined for 


corresponding elements of vector double x. A vector double is returned that contains the value of ex, for x,, such 
that: 


e x, == frac* FLT RADIX™ 
° | frac] is in the interval [1, FLT RADIX) 
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If x; is O, the corresponding result is undefined and no error is reported. 





On the SPU, if x; is O, the corresponding result is -HUGE VALF. If x; is infinite, the corresponding result is postive 
infinite. If x; is a NaN, the corresponding result is also a NaN. 


logbf4: Represent Float as Fraction Greater Than 1 and Exponent 


(vector float) logbf4 (vector float x); 


An integer exponent ex and a fraction frac, that represents the value of a finite element are determined for 
corresponding elements of vector float x. A vector float is returned that contains the values of ex, for x;, such that: 


e x;-- frac * FLT RADIX^ 
° | frac] is in the interval [1, FLT RADIX) 


If x; is O, the corresponding result is undefined and no error is reported. 


modfd2: Represent Double as Proper Fraction and Exponent (SPU Only) 


(vector double) modfd2 (vector double x, vector double *pint); 


Each element of vector double x is split into an integral part i, and a fractional part frac. A vector double is 
returned that contains the corresponding frac, elements, and another vector double is stored in *pint that contains 
the corresponding i, elements, such that: 


. x; == fracti 
° | £rac| is in the interval [0, 1) 
e both frac and i have the same sign as x, 


modff4: Represent Float as Proper Fraction and Exponent 


(vector float) modff4 (vector float x, vector float *pint); 


Each element of vector float x is split into an integral part i, and a fractional part £rac;. A vector float is returned 
that contains the corresponding frac; elements, and another vector float is stored in *pint that contains the 
corresponding i; elements, such that: 


. xi == frac + i 
° | £rac| is in the interval [0, 1) 
e both frac and i have the same sign as x; 
nearbyintd2: Find Nearest Integer for Double (SPU Only) 
(vector double) nearbyintd2 (vector double x) 
A vector double is returned that contains the corresponding elements of vector double x rounded to the nearest 


integer consistent with the current rounding mode, but without raising an inexact floating-point exception. 


nearbyintf4: Find Nearest Integer for Float 


(vector float) nearbyintf4 (vector float x) 


A vector float is returned that contains the corresponding elements of vector float x rounded to the nearest integer, 
consistent with the current rounding mode, but without raising an inexact floating-point exception. 


On the SPU, the rounding mode for a float is always towards zero. 


negated2: Negate Double (SPU Only) 
(vector double) negated2 (vector double x); 





A vector double is returned that contains —x, for corresponding elements of vector double x. 
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negatef4: Negate Float 


(vector float) negatef4 (vector float x); 


A vector float is returned that contains —x; for corresponding elements of vector float x. 


negatei4: Negate Signed Integer 
(vector signed int) negatei4 (vector signed int x); 
A vector signed int is returned that contains —x, for corresponding elements of vector signed int x. 


If —x, cannot be represented, the corresponding result is undefined and no error is reported. 


negatell2: Negate Signed Long Long Integer (SPU Only) 


(vector signed long long) negatell2 (vector signed long long x); 
A vector signed long long is returned that contains —x; for corresponding elements of vector signed long long x. 


If —x; cannot be represented, the corresponding result is undefined and no error is reported. 


nextafterd2: Find Next Integer After for Double (SPU Only) 





(vector double) nextafterd2 (vector double x, vector double y) 


A vector double is returned that contains the next representable value after x; in the direction of y, for corresponding 
elements of vector double x and vector double y. If x; is equal to y, the result is y;. 


If the magnitude x; is the largest finite value representable, the result is undefined. 


nextafterf4: Find Next Integer After for Float 


(vector float) nextafterf4 (vector float x, vector float y) 





A vector float is returned that contains the next representable value after x, in the direction of y, for corresponding 
elements of vector float x and vector float y. If the element of x; is equal to y,, the result is y;. 


If the magnitude of x; is the largest finite value representable, the result is undefined. 


powd2: Raise Double to Double Power (SPU Only) 


(vector double) powd2 (vector double x, vector double y); 


A vector double is returned that contains x, raised to the power of y;, x”,, for corresponding elements of vector 
double x and vector double y. 


If x, is finite and negative and y, is finite and not a integer value, the corresponding result is undefined and no error is 


reported. 


powf4: Raise Float to Float Power 


(vector float) powf4 (vector float x, vector float y); 


A vector float is returned that contains x, raised to the power of y,, x";, for corresponding elements of vector float x 
and vector float y. 


On the SPU, if the result would be greater than HUGE VALF, the result is saturated to HUGE VALF and no error is 
reported. 
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recipd2: Reciprocal of Double (SPU Only) 


(vector double) recipd2 (vector double x); 
A vector double is returned that contains the reciprocal of the corresponding elements of vector double x. 
The function handles special cases as follows: 


e When x, is «Inf, the result is 0 with the sign of x; 
e When x; is 0, the result is Inf with the sign of x; 


e When x; is NaN, the result is NaN 
recipf4: Reciprocal of Float 
(vector float) recipf4 (vector float x); 
A vector float is returned that contains the reciprocal of the corresponding elements of vector float x. 
The function handles special cases as follows: 


e When x; is «Inf, the result is 0 with the sign of x; 
e When x, is 0, the result is HUGE VALF with the sign of x; 





e When an element of x, is NaN, the result is NaN 


remainderd2: Remainder of Doubles (SPU Only) 





(vector double) remainderd2 (vector double x, vector double y); 


A vector double is returned that contains the remainder x; REM y; for the corresponding elements of vector double x 
and vector double y. 





If y; is zero, the corresponding element of the result is undefined and no error is reported. 


remainderf4: Remainder of Floats 


(vector float) remainderf4 (vector float x, vector float y); 








A vector float is returned that contains the remainder x; REM y;, for the corresponding elements of vector float x and 
vector float y. 


If y; is zero, the corresponding element of the result is undefined and no error is reported. 


remquod2: Remainder Function of Double (SPU Only) 
(vector double) remquod2 (vector double x, vector double y, vector signed long long 
*pquo) 
This function returns the same vector double result as remaindera2 (). In addition a vector signed long long is 
stored in *pquo that contains the corresponding element values whose sign is the sign of x; / y; and whose 
magnitude is congruent modulo 2^ to the magnitude of the integral quotient of x; / y;, where n is an implementation- 
defined integer greater than or equal to 3. 


remquof4: Remainder Function of Float 


(vector float) remquof4 (vector float x, vector float y, vector signed int *pquo) 

This function returns the same vector float result as remainderf£4 (). In addition a vector signed int is stored in 
*pquo that contains the corresponding element values whose sign is the sign of x; / y, and whose magnitude is 
congruent modulo 2” to the magnitude of the integral quotient of x; / y,, where n is an implementation-defined integer 
greater than or equal to 3. 
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rintd2: Round Double to the Nearest Integer (SPU Only) 


(vector double) rintd2 (vector double x); 


A vector double is returned that contains the corresponding elements of vector double x rounded to the nearest 
integer, consistent with the current rounding mode. 


rintf4: Round Float to the Nearest Integer 


(vector float) rintf4 (vector float x); 


A vector float is returned that contains the corresponding elements of vector float x rounded to the nearest integer, 
consistent with the current rounding model. 


On the SPU, the rounding mode for float is always towards zero. 


roundd2: Round Double (SPU Only) 


(vector double) roundd2 (vector double x); 


A vector double is returned that contains the rounded elements of vector double x. Rounding is done to the nearest 
integer value in floating-point format. Halfway cases are rounded away from zero regardless of the current rounding 


direction. 
roundf4: Round Float 


(vector float) roundf4 (vector float x); 


A vector float is returned that contains the rounded elements of vector float x. Rounding is done to the nearest integer 
value in floating-point format. Halfway cases are rounded away from zero regardless of the current rounding direction. 


rsqrtd2: Reciprocal Square Root of Double (SPU Only) 


(vector double) rsqrtd2 (vector double x); 


A vector double is returned that contains the reciprocal of the square root of x, for the corresponding elements of 
vector double x. Special cases are handled as follows: 


e When x, is less than 0, the result is NaN 

e When x; is «Inf, the result is +0 

e When x; is 0, the result is Inf with the sign of x; 
e When x; is NaN, the result is NaN 


rsqrtf4: Reciprocal Square Root of Float 


(vector float) rsqrtf4 (vector float x); 


A vector float is returned that contains the reciprocal of the square root of x; for the corresponding elements of vector 
float x. Special cases are handled as follows: 


e When x, is less than 0, the result is NaN 

e When x; is +Inf, the result is +0 

e When x, is 0, the result is Inf with the sign of x, 
e When x, is NaN, the result is NaN 


On the SPU, if x, is less than zero, the corresponding result is undefined. 


scalblind2: Scale Double by Long Long Integer (SPU Only) 


(vector double) scalbllnd2 (vector double x, vector signed long long n) 


A vector double is returned that contains x, efficiently multiplied by 2^ for corresponding elements of vector double x 
and vector signed long long n. 
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scalbnf4: Scale Float by Integer 
(vector float) scalbnf4 (vector float x, vector signed int n) 


A vector float is returned that contains x; efficiently multiplied by 2” for corresponding elements of vector float x and 
vector signed int n. 


signbitd2: Sign Bit of Double (SPU Only) 


(vector unsigned long long) signbitd2 (vector double x) 


A vector unsigned long long is returned that contains the elements defined below for corresponding elements of 
vector double x. 


e All bits of the resulting element are set to 1 if the sign bit is set in x; 
e Zero otherwise 


signbitf4: Sign Bit of Float 
(vector unsigned int) signbitf4 (vector float x) 


A vector unsigned int is returned that contains the elements defined below for corresponding elements of vector 
float x. 


e All bits of the resulting element are set to 1 if the sign bit is set in x, 
e Zero otherwise 


sincosd2: Sine and Cosine of Double (SPU Only) 


(void) sincosd2 (vector double x, vector double *sx, vector double *cx); 


A vector double is stored in *sx and a vector double is stored in *cx that contain the respective sines and cosines of 
the corresponding elements of vector double x. 


The results of sincosd2 () may not be accurate for very large values of x;, and no error is reported. 
Implementations should document the point at which accuracy is lost. 


sincosf4: Sine and Cosine of Float 
(void) sincosf4 (vector float x, vector float *sx, vector float *cx); 


A vector float is stored in *sx, and a vector float is stored in *cx, that contains the respective sines and cosines of 
the corresponding elements of vector float x. 


The results of sincos£4 () may not be accurate for very large values of x,, and no error is reported. 
Implementations should document the point at which accuracy is lost. 


sind2: Sine of Double (SPU Only) 
(vector double) sind2 (vector double x); 
A vector double is returned that contains the corresponding sines of the elements of vector double x. 


The results of sind2 () may not be accurate for very large values of x,, and no error is reported. Implementations 
should document the point at which accuracy is lost. 


sinf4: Sine of Float 
(vector float) sinf4 (vector float x); 
A vector float is returned that contains the corresponding sines of the elements of vector float x. 


The results of sin£4() may not be accurate for very large values of x;, and no error is reported. Implementations 
should document the point at which accuracy is lost. 
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sinhd2: Hyperbolic Sine of Double (SPU Only) 


(vector double) sinhd2 (vector double x) 


A vector double is returned that contains the corresponding hyperbolic sines of the elements of vector double x. 


sinhf4: Hyperbolic Sine of Float 


(vector float) sinhf4 (vector float x) 


A vector float is returned that contains the corresponding hyperbolic sines of the elements of vector float x. 








On the SPU, element values of the result that are greater than HUGE VALF are returned as HUGE VALF, and no error 
is reported. 


sqrtd2: Square Root of Double (SPU Only) 
(vector double) sqrtd2 (vector double x); 


A vector double is returned that contains the real square roots x^/?; for the corresponding elements of vector 
double x. 


This function handles special cases as follows: 


e When x; is less than 0, the result is NaN 
e When xi is +Inf, the result is +Inf 
e When x, is 0, the result is 0 with the sign of x; 


e When x, is NaN, the result is NaN 
sqrtf4: Square Root of Float 
(vector float) sqrtf4 (vector float x); 
A vector float is returned that contains the real square roots x^? for the corresponding elements of vector float x. 
This function handles special cases as follows: 


e When x; is less than 0, the result is NaN 

e When x, is +Inf, the result is +Inf 

e When x; is 0, the result is 0 with the sign of x; 
e When x, is NaN, the result is NaN 


On the SPU, the result is undefined when x; is negative. 


tand2: Tangent of Double (SPU Only) 
(vector double) tand2 (vector double x); 
A vector double is returned containing the corresponding tangents of the elements of vector double x. 
The results may not be accurate for very large values of x, and no error is reported. Implementations should 


document the point at which accuracy is lost. 


tanf4: Tangent of Float 


(vector float) tanf4 (vector float x); 
A vector float is returned containing the corresponding tangents of the elements of vector float x. 


The results may not be accurate for very large values of x, and no error is reported. Implementations should 
document the point at which accuracy is lost. 
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tanhd2: Hyperbolic Tangent of Double (SPU Only) 


(vector double) tanhd2 (vector double x) 


A vector double is returned that contains the corresponding hyperbolic tangents of the elements of vector double x. 


tanhf4: Hyperbolic Tangent of Float 
(vector float) tanhf4 (vector float x) 


A vector float is returned that contains the corresponding hyperbolic tangents of the elements of vector float x. 


tgammad2: Gamma of Double (SPU Only) 


(vector double) tgammad2 (vector double x) 


A vector double is returned that contains the corresponding results of the gamma function applied to the respective 
elements of vector double x. 


If x; is a negative integer, the corresponding element of the result is undefined and no error is reported. 


tgammaf4: Gamma of Float 


(vector float) tgammaf4 (vector float x) 


A vector float is returned that contains the corresponding results of the gamma function applied to the respective 
elements of vector float x. 


If x; is a negative integer, the corresponding element of the result is undefined and no error is reported. 


truncd2: Truncate Double (SPU Only) 


(vector double) truncd2 (vector double x); 


A vector double is returned that contains x; rounded to the nearest integer n that is not larger in magnitude than x; 
(rounded towards zero) for each corresponding element of vector double x. 


truncf4: Truncate Float 


(vector float) truncf4 (vector float x); 


A vector float is returned that contains x, rounded to the nearest integer n that is not larger in magnitude than x, 
(rounded towards zero) for each corresponding element of vector float x . 
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